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1 Introduction 



Quark and lepton masses and mixing angles are free parameters of the Standard Model 
(SM). They are known with various degree of accuracy, from the precision direct mea- 
surements of charged lepton masses, a combination of experimental data and theoretical 
arguments for quark masses and relatively precise measurements of quark mixing angles 
and the phase parameterizing the violation of CP invariance in the SM. A new element in 
this picture are neutrino masses and mixing angles. There is at present strong experimen- 
tal evidence for neutrino oscillations whose most obvious and most natural explanation 
is that neutrinos have non-zero masses and the neutrino mass eigenstates are different 
from the weak interaction eigenstates. Although far from having final interpretation, 
the present experimental data give interesting preliminary information about the neu- 
trino mass sector and the forthcoming experiments are expected to resolve the remaining 
ambiguities. 

On the theoretical side, the origin of the interactions giving rise to fermion masses 
is a problem that cannot be addressed in the framework of the SM. The physical scale 
(we shall call it M) of the, still unknown, theory of fermion masses is certainly above the 
electroweak scale and quite likely it may even be close to the GUT or Planck scales. Such 
a high scale is suggested, for instance, by the see-saw interpretation of the magnitude of 
the neutrino masses indicated by experiment. Ultimately, the theory will predict at the 
scale M the running fermion mass parameters (and perhaps also other parameters) of the 
effective low energy theory, describing physics at energies < M. [| 

To relate the mass parameters of the effective low energy theory to the experimentally 
measured quantities one has to include quantum corrections that already do not depend on 
the specific theory of fermion masses but only on the low energy effective theory. Clearly, 
close to the electroweak scale this theory is the Standard Model (SM). One possibility is 
that the SM remains the correct effective theory up to the scale M. This is conceivable 
particularly if the scale M is relatively low. Another possibility is that the SM needs to be 
embedded into a bigger effective theory already much below M. The latter case may, for 
instance, happen if the low energy supersymmetry is realized in Nature. In this review we 
discuss quantum corrections to neutrino masses both in the SM and in its supersymmetric 
extension, the Minimal Supersymmetric Standard Model (MSSM). 

Quite generally one can distinguish two classes of quantum corrections that enter on 
a somewhat different footing. The first one is given by the RGEs describing the evolution 

^Thc running mass parameters of the effective low energy tlieory can in principle be calculated in 
the underlying theory at any renormalization scale Q by proper inclusion of the high energy threshold 
corrections. Choosing Q k M one minimizes those corrections. 
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of the fermion mass parameters in the SM or MSSM from the scale M down to some 
scale close to the electroweak scale. Another source of quantum corrections are the so- 
called low energy threshold effects. Strictly speaking these are the corrections necessary 
to express measurable quantities like neutrino masses and mixing angles in terms of the 
running (renormalized) parameters of the effective theory Lagrangian. Formally, their 
inclusion renders the prediction for observables independent of the choice of the scale to 
which the RGEs are integrated. In the SM the threshold corrections are unambiguous 
and, if that scale is taken to be close to M^, they are small and can be neglected. Such 
corrections may, however, be very important in the MSSM since there they depend also 
on the sparticle masses and couplings. 

In Section 2 we review the neutrino masses and mixing and stress the differences with 
the quark sector. Contrary to the small mixing in the latter one, at least one neutrino 
mixing angle and, quite likely two, are close to maximal A maximal or bimaximal 
mixing would be quite natural for (approximately) degenerate masses, much larger than 
their differences 0. Such mass patterns are very different from the hierarchical masses 
of the charged fermions but are consistent with experiment. Indeed, only mass squared 
differences can be inferred from the data and not the neutrino masses themselves. 

The observed (mixing) and potential (mass pattern) differences between the neutrino 
and quark sectors provide strong motivation for studying quantum corrections to neu- 
trino masses and mixing. It is well known that the corrections are small for the quarks, 
just because of their hierarchical masses and mixing. As we shall review in this paper, 
that remains true for neutrinos if their masses are hierarchical too. However, quantum 
corrections may give strong, qualitatively new, effects if the neutrino masses are not hi- 
erarchical. 

The main parts of this review are Sections 4 and 5. In Section 4 we discuss quantum 
corrections described by the renormalization group evolution and in Section 5 the potential 
effects of the low energy threshold corrections are reviewed. 

A brief overview of quantum corrections in the neutrino sector is given in Section 6. 
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2 Fermion masses and mixing 



For quarks and charged leptons the Particle Data Group gives the following values of 
the masses^: 

m„ = 1.5-5.0 MeV, = 1.15 - 1.35 GeV, mj = 174.3 ± 5.1 GeV, 

ma = 3.0 - 9.0 MeV, = 60 - 170 MeV, = 4.0 - 4.4 GeV, (2.1) 

me = 511 keV, = 105.7 MeV, = 1.777 GeV 

and of the quark mixing (absolute values of the elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix): 

/0.9742 - 0.9757 0.219- 0.226 0.002 - 0.005 \ 

0.219 - 0.225 0.9743- 0.9749 0.037- 0.043 (2.2) 

^ 0.004- 0.014 0.035 - 0.043 0.9990 - 0.9993 y 
The pattern of quark masses and mixing is clear: hierarchical masses and small (and 
also) hierarchical mixing: the larger the mass difference, the smaller the mixing angle. In 
the framework of the SM, fermion masses and mixing angles are used to fix the Yukawa 
matrices which are the actual free parameters of the theory. They cannot, however, 
be uniquely reconstructed from the experimental data: there is always the ambiguity of 
rotating the electroweak basis and all we can have at present is merely a phenomenological 
parametrization. 

Neutrino masses and mixing are inferred from a number of oscillation experiments. 
The interpretation of neutrino oscillations in terms of massive neutrinos is the most nat- 
ural one. To discuss this interpretation, let us first recall the oscillation pattern for two 
hypothetical neutrinos with masses mi and 1712 whose quantum fields are linear combina- 
tions of two quantum fields ua and ub that are weak eigenstates 

z/i = cos {^ua + sin 'duB 

z/2 = —sin'dh'A + cos{}h'B (2.3) 

The probabilities that in a given experiment a neutrino produced in the interaction of 
the charged lepton of fiavour A with the boson is detected at a distance L as the 
neutrino creating the charged lepton of fiavour B, B A, or A are0 

o2x2/ N •2oQ-2 / -f'cxp 

B [ua —>■ Vb] ~ ^-'J— ^ ' 




P'^'^yA - ^a) = 1 - P'^'^yA - ys) (2.4) 



^For M, d and s quarks the quoted values are the running masses in the MS scheme at Q = 2 GeV. 
The values given for the c and h quarks are the MS running masses at Q = TOc and Q = mb, respectively. 
Finally mt given here is the pole mass. 

^The formulae like (2^) are usually derived in the framework of quantum mechanics. The proper 
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We see that the oscillation probability depends on Am^ = m^^ — m^^, the distance Lexp, 
the mixing angle "i? and the neutrino energy E^xp- It is clear that for observing the 
oscillation pattern the factor Am'^L/4E should be of order 0{1). Q Hence, in general, 
longer oscillation distances are necessary to probe smaller mass squared differences. On 
the other hand, for Am'^L/AE ^ 1 the transition probability is P(i^a ^b) ~ |sin^2^9, 
i.e. it is insensitive to Am^, because in realistic applications the expression (|2.4|) has to 
be averaged over some non-zero interval of the initial neutrino energies AE 

We can now summarize the results of various experiments. They are usually inter- 



preted in terms of the effective 2x2 parametrization ( |2.4| ). Historically, the first infor- 
mation came from the experiments measuring the flux of neutrinos produced in the 
Sun (for references see e.g. the review 0). A convincing evidence for the solar neutrino 
oscillation was provided by the Kamiokande and Superkamiokande experiments which es- 
tablished a strong 50%) suppression of the flux of neutrinos from the nuclear reaction 
^Be* + e+ + z/g. The most plausible explanation of those results is the transmu- 
tation of i/g's produced in the core of the Sun into another type of neutrinos such as z/^, 
Ut- and/or the so-called sterile neutrino t'steriie which does not interact with the or 
bosons. Such a transmutation can occur either during their flight from the Sun to the 
Earth (the so-called vacuum oscillations (VO)) or through the resonant transition in the 
matter of the outer layers of the Sun (the so-called MSW effect [^). As follows from the 
formulae ( p.4| ) with A = e and 5 = /x, r or s, in the case of VO with L 1.5 x 10® km (the 
Sun-Earth mean distance) and for mean ^B neutrino energy of order ~ 10 MeV, the 
solar neutrino experiments are sensitive to Aml^i ~ O{10~^^^~^^^) eV^. The deficit is then 
explained with sin^ 2t9soi > 0.7. In the case of the MSW resonant conversion the formulae 
( p.4|) for the transition and survival probabilities are replaced by more complicated ex- 
pressions which depend also on the electron and neutron number densities in the Sun. 
The observed neutrino deficit can be then explained either for Am^^j ~ O{10~^) eV^ 
and sin2 279^01 ~ C(10-(3-2)) (the so-called SAMWS solution) or for Am^^ ~ 0(10"^) eV^ 
and sin^ 279801 > 0.5 (the so-called LAMWS solution). 



picture of neutrino oscillations is however the field theoretical one Q|. Due to the field mixing (2.3) 
there is in general a nonzero amplitude for emission (absorption) of any of the neutrino mass eigenstates 

in the interaction of charged lepton of flavour A with the VF-boson. The change of the neutrino 
flavour (inferred from the flavours of the charged leptons) is in this picture due to the coherent sum of 
Feynman diagrams describing exchanges of all virtual mass eigenstates of neutrinos between the emission 
and detection vertices. Taking properly into account the effects of wave packets describing initial and 
final states it can be shown that the formula (2.4) is in most cases a sufficient approximation to the 
full result, which automatically accounts for decoherence effects and depends on the overlap of the wave 
packets describing the initial and final states 

^Recall that ^ = 1.27 x (^"'{^J^^evf 
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The deficit of neutrinos was also revealed by the measurements of the flux of 
and z/^ neutrinos produced together with i/g's in the Earth's atmosphere by the cosmic 
rays. The results of the Superkamiokande experiment [0 are most easily explained by 
oscillation into another type of neutrinos. The oscillatory explanation is further supported 
by the zenith angle dependence of the z/^ and z/^ flux deficit. For typical z/^ energies of 
order ~GeV and 20 km < L < 1.3 x 10^ km, the observed deficit is explained for 
AMftj^ ^ 3.2 X 10~^ eV^ and sin^ 2'(9atm > 0.82. Also, the Superkamiokande data seem to 
favour z/^ — > u^- oscillations over — > z/gteriic |@]- 

Neutrino oscillations are also intensively searched for in various reactor or accelerator 
based experiments (for review see R P]). Except for the LSND experiment reporting [|IU| a 



positive signal for the z/^ Ue oscillations, which would require 0.1 eV^ < AM^ < 1 eV^ for 
the oscillatory explanation, the results of all other experiments are negative. Particularly 
strong constraint comes from the CHOOZ reactor experiment. It excludes disappearance 
of z/g with mean energies ~ 3 MeV at the distance L ^ 1 km. In the 2x2 oscillation 
framework this translates into the limit 

sin2 279,eact < 0.1 for AM^ > 10'^ eV^ (2.5) 

Except for the unconfirmed LSND result, all experiments are consistent with oscil- 
lations between the three known neutrino flavours. Moreover, recent SNO measurement 
pj| of the total solar neutrino flux (i.e. the flux of z/g, z/^ and Vr) combined with the 



Superkamiokande data |T2[ strongly supports this assumption. Therefore, in the rest of 
the review we assume that the neutrino sector of the low energy effective theory consists 
of 3 active neutrinos only. 

Within the 3x3 framework the solar, atmospheric and reactor data become interrelated 
and one can draw more definite conclusions. To make the discussion of the neutrino 
experiments more complete we recall the formula for the transition probabilities P{va 
Ub) for the case of 3 neutrinos: 



m, , L 



Aa 



= E lUnaUXaf + Re E E f/Baf/^f/l.f/A.e^^™-^/'^ (2.6) 

where Am^f, = mf,^ ~''^ut- "^^^ complex matrix elements are defined in the basis in which 
charged lepton mass matrix is diagonal, by the decomposition of the neutrino fields 
belonging to SUl{2) doublets into the mass eigenstates field: 

l^A = T.UAaiya. (2.7) 
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For a real matrix U the expression ( p^.6D can be written as 



P{ua ub, L, E) = 6ab-4: J2 UAaUBaUAbUsb sm^ ( 

a>b=l \ 



AE 



To simplify the notation we will denote: Am = Amg^j = \ra2~m{\, AM = AM^^^ = 
Img — mil- For Am^ ^ AM^ the measurements of the experiments listed above can now 
be summarized as follows: 



Am^Lsoi 



sol 



/ AM^Latm 



PatmK-^M) = l-4(l-[/4)[/4sin^(^ 4eJ j ^^-^^ 



^'reactli^e ^ Z^e) = 1 - 4 (l - f/^g) ^J^g siu^ 



fAM^L 



react 



A convenient parametrization for the mixing matrix U is (see also the next Section): 



MNS 



C12C13 S12C13 Si^e 

-S12C23 - Ci2S23Si3e*'' C12C23 - Si2S23Sl3e*'' S23C13 
V ■S12S23 — Ci2C23Si3e*'^ — C12S23 — Si2C23Si3e*'^ C23C13 / 



X $ (2.10) 



where Cij = cos 6ij {sij = sin 6ij) and $ is a diagonal matrix $ = diag{e^°'^,e''"^, 1) (or a 
similar diagonal matrix with two phases and $11 = 1 or $22 = !)• Our convention is such 
that for vanishing mixing angles dij we have ui = Ug, 1/2 = and u-^ = 1/^, i.e. no ordering 
of the neutrino masses rrii,^ is assumed. CP is conserved if 5 = mod vr and ai^2 = mod 
7t/2. The relation to the effective 2x2 parametrization then reads 

sin^ 2t9i.eact = sin^ 26*13 

2 n 2/1/2/1 sin^ 26^13 \ 

sm 2'ds,tm = sm 26'23 cos fc'13 + — ^ 

\ 4cosi3 / 

sin^ 2?9soi ~ sin^ 2^12 (l - 2 sin^ ^13) (2.11) 



In the limit of small ^13 angle imposed by the CHOOZ result ( p.5| ) (non-negligible entries 
U23 and ?733 and unitarity of U exclude the possibility of |6'i3| ^ 7r/2) we get 

sin^ 2?9atm ~ sin^ 26*23 

sin2 2^?soi ~ sin2 2^i2 (2.12) 

Finally we mention that the neutrino Majorana mass term of the form £niass = 
— ^rn^^UA^B + H.c violates the lepton number conservation and leads to the neutrino- less 
double beta decay. Non-observation of such decays implies |TBI 



< 0.35(0.27) eV at 90% (68%) C.L., (2.13) 
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i.e. 

|m^,c?2C?3e^'"' +^^2S?2Ci3e^*"' +m^3S^3e2*'^| < 0.35(0.27) eV at 90% (68%) C.L.(2.14) 

Concluding this section, we stress that the pattern of neutrino mixing is distinctly 
different from the one observed in the quark sector. At least one mixing angle {623) is 
large and the recent data favour solutions with two large mixing angles The data 
indicate Am^ <^ AM^ but, since the masses are not measured, several options for the 
values of the masses themselves are still possible: 

i) AM^ ^ ml^ > "^^2(1) ^ "^^1(2) (hierarchical) 

ii) AM^ ^ ml^^ml^ ml^ or AM^ ^ ml^ ■ml_^,ml^, with m1_^^m1^ ^ Am^ (partly 
degenerate) 

in) ~ ^ m^^ and all of order or larger than AM^ (degenerate). 

Special cases of the partly degenerate and degenerate patterns are equal masses, jm^J = 
Im^J or \my^ \ = \m^^\ = \mi,^\, respectively. In our language we refer to them as two-fold 
and three-fold degeneracies, to distinguish them from, more general, partly degenerate 
and degenerate pattern defined by ii) and Hi). The last two possibilities are very dif- 
ferent from the pattern known from the quark sector but almost maximal (bimaximal) 
mixing suggested by the experimental data makes them an interesting alternative to the 
hierarchical pattern. It is important to remember that in that respect neutrinos may 
be qualitatively different from quarks because, unlike the charged fermions, they can be 
Majorana particles. As we shall see, the magnitude and the importance of quantum 
corrections depends strongly on the assumed pattern of the masses. 



3 Neutrino masses in the effective theory 

We begin by recalling how the quark masses are incorporated in the SM (or the MSSM) 
Lagrangian. We shall use the Weyl spinor notation (see e.g. |0), which is particularly 
convenient for Majorana particles. The relation to the standard Dirac notation is ex- 
plained in the Appendix A. The original Lorentz and SU{3) x SU{2) x U{1) invariant 
Lagrangian contains the Yukawa interactions of the form 

£Yuk = -e,,H,u'^^Y^^qf - H^d^^Vj^f - H^e^^Y^^lf + H.c. (3.1) 

where Yu, Yd, Y e are a priori arbitrary complex Yukawa matrices, ei2 = —£21 = —1, the 
upper case letters enumerate the three generations of the matter fermions 

(T, e\ and / (3.2) 
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described by the left-handed Weyl spinors and H is the SM Higgs doublet. In the MSSM 
there are two Higgs doublets, if*-"^ and and one has to replace: h\^'^ and 

H* eijH\'^\ Replacing the Higgs doublet H (doublets H^""^ and H^"^^ in the MSSM) 
by its vacuum expectation value leads to the following mass terms in the SM (MSSM) 
Lagrangian: 

Cr^^s = -u''''{vY^)^%^ - d^^{vY,)^''d^ - e'^{vY,)^''e^ + H.c. (3.3) 

(In the MSSM vY^ v^Y^, vY ~^ ~VdY fi,e-) The fermion mass matrices are diago- 
nalized by the unitary chiral rotations 

d DLd, d^ d'D^ji, (3.4) 
e ^ Elc, ^ e'^E^j,, 

which give: 

U^j^Yy^UL = diag(?/„, y^, yt) 

D^RYdDL = diag(?/d, y^, yt) (3.5) 
E^j^Y^El = diag(ye, 2/^' Vr)- 

The only remnant of the nontrivial Yukawa matrices is then the CKM matrix 

V^CKM = V = U{Dl (3.6) 

appearing in the interactions of the charged quark currents with (massive) charged vector 
bosons W^. It can be shown (see e.g. [^) that the CKM matrix can be parameterized 
by 3 angles and one phase. The other 6 — 1 = 5 phases can be absorbed in the redefinition 
of the quark fields. 

In the SM defined as a renormalizable theory, with the neutrino and Higgs fields 
transforming as components of doublets of the SUl{2) gauge symmetry, neutrinos remain 
massless. Moreover, as a consequence of the gauge symmetry, renormalizability and the 
field content, the Lagrangian has two accidental global U{1) symmetries that ensure 
baryon and lepton number conservation. In the absence of any mass terms for the neutrino 
fields z/^, one can always perform the rotation z/ V^z/ with Vl = El so that the 
counterpart of the CKM matrix in the leptonic sector is trivial. 

There are two easy possibilities for extending the particle content of the SM so that 
neutrino mass is generated by renormalizable interactions. One is to couple two lepton 
SUl{2) doublets to a SUl{'2.) Higgs triplet (singlets would violate the electric charge 
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conservation). This possibility is not particularly attractive for several reasons. The 
smallness of the neutrino masses would require either ftripiet ^ "^doublet or the triplet 
couplings orders of magnitude smaller than the other Yukawa couplings. Moreover, the 
introduction of the triplet would make the parameter p a free parameter of the theory 
(hence not calculable) with its own counterterm allowing to adjust its value at will. The 
other possibility is to introduce a number of new SUl{2) x f/y(l) singlet left-handed 
leptonic fields z/^ to the Lagrangian and couple them to the Higgs doublet in the same 
way as the quark fields u^: 



where K = 1, . . . , K^ax and in principle Kmax could be arbitrary. The interaction (|3.7|) 
preserves the known structure of the SM and, if = 1, 2, 3, the symmetry between quarks 
and leptons is restored. For singlet fields v'^ the Majorana mass term 



can also be added to the Lagrangian. In general both terms should be included in the the- 
ory defined by the SUl{2) x f/y(l) gauge symmetry and renormalizability. The Majorana 
mass term can be, however, eliminated by imposing the additional global U{1) symmetry 
ensuring conservation of the lepton number L. The fields v'^ must then have L = —1, 
opposite to L of the leptonic doublets k. Their complex conjugate u'^ (see Appendix A) 
can be then interpreted as the right-handed neutrinos. Neutrinos are then Dirac particles 
like the other fermions. There are two reasons why Dirac masses are not so attractive. 
One is, again, the need for very small numerical values of the Yukawa couplings Y^. The 
other is that the lepton number conservation has to be imposed as an additional global 
symmetry (remember that in the SM it is a consequence of the field content and of the 
renormalizability and not an additional assumption). 

The presence of the Majorana mass term, which breaks the global U{1) symmetry, 
inevitably makes neutrinos Majorana particles (in fact what allows to interpret in the SM 
the two helicity states described by the u field as a particle and an anti-particle is just the 
lepton number!) Taking for simplicity threeQ singlet neutrinos and replacing H {H^'^^ 
in the MSSM) by its vacuum expectation value, the general form of the neutrino mass 
matrix (in the absence of Higgs triplets) is 



(3.7) 




(3.8) 




(3.9) 



■'The formula (3.9) is easily generalized to an arbitrary number of singlet neutrinos j/"^. 
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where = vY,y and MMaj are 3x3 matrices in the basis {u,!/'^). If both these ma- 
trices are diagonal with eigenvalues and Ma, respectively, the physical neutrinos are 
mixtures of ua and z/^ with the masses 

1 



m,A± = - (^Ma ± sjMl + Am\j . (3. 10) 

In general, each of the six leptonic fields ua and z/^ is a linear combination of the six 
Majorana neutrinos. There is the possibility of oscillations u —>■ v'^. The v'^ neutrinos 
are singlets of the SUl{2) x f/(l) gauge group and are sterile neutrinos.^ The situation 
simplifies for M ^ v where M is the overall scale of the MMaj entries. For simplicity, we 
identify it with the introduced earlier scale M, although they could actually be different 
scales. In this limit there are two sets of mass eigenvalues. They are of the order 0{M) 
and C(f^/M), respectively, and three of the mass eigenstates are very heavy and decouple 
from the physics at the electroweak scale. This is the so-called see-saw mechanism [|l^ 
[] The oscillations occur then effectively among the three "active" neutrinos u (the i/'^'s 
composition is dominated by heavy eigenstates). 

It is instructive to repeat the decoupling procedure in the field theoretical language. 
The effects of the Born diagram. Fig. present in the full theory valid above the scale 
M are, up to C(l/M^), reproduced in the effective theory, describing physics below the 
scale M, by a non-renormalizable operator of dimension 5 [T^|: 



A£ = --;^C^''{e.,Hdf){eiMj^) + H.c. (3.11) 



4M 
in which 

M-'C = Yl{Mu.j)-'Y,. (3.12) 

The overall scale M of the entries of the matrix Muaj has been factorized out from C to 
make it dimensionless. Thus, if the scale M is high enough, the effects of the Majorana 
masses can be discussed in the SM or MSSM supplemented by the operator ( p.ll| ). We 
should also stress that this operator is the only one of dimension 5 contributing to the 
neutrino Majorana mass matrix. Other possible contributions are of higher dimension 
pO| . Thus, one may expect that in the effective SM or MSSM the neutrino masses are 
indeed described by the operator ( |3.11|) , even if its origin is different from the see-saw 



^Note that such oscihations exist only for Majorana neutrinos and are forbidden by L conservation for 
Dirac neutrinos i.e. the right-handed neutrinos P"^ have correct L but wrong chirahty to osciUate into;^. 

^However in some cases, one of the formally 0{M) mass eigenvalues can remain small as a result of 
some symmetries of the underlying theory of neutrino masses. The corresponding light state can play 
the role of a sterile neutrino 
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mechanism. In particular, degenerate or partly degenerate patterns of neutrino masses 
require some interplay between the parameters entering the formula ( p.l2| ), but one can 
also imagine that such patterns originate from other mechanisms at the scale M. 




Figure 1: Diagram generating the dimension 5 operator (|3.11|) . 

We end this section with the discussion of the parametrization of the neutrino mixing 
matrix. In the effective theory, after the electroweak symmetry breaking, the operator 
is the source of the Majorana mass for the three active neutrinos z/, £mass = 
v^v^ + H.c. This mass matrix can be diagonalized by an additional unitary 
rotation [/^""Ua- Recall we work in the basis in which the leptonic Yukawa coupling 

is already diagonal and the matrix C is assumed to be given in the same basis. We have: 

— (U^CUY' = —C-S^' = 5^' (3.13) 

The matrix U will therefore appear in the couplings of neutrinos to bosons and 
charged Goldstone (and Higgs) boson.0 Being unitary, the matrix U depends on 3 angles 
and 6 phases and can be conveniently written as U^"' = e^'^'^f/j^^s where Umtsss is given 
in eq. (|2.10|) . Contrary to the quark case, if the Majorana masses m^^ in eq. ( |3.13|) 
are to be real and positive, the only freedom that remains is the possibility to re-phase 
independently the three Dirac fields of the charged leptons ipeA ~^ Gxp{iipA)i'eA 
eA — exp(zy9A)eA for the left-handed and — ^ exp(— 2(/9A)e$i for their right-handed 
components). It is then the MNS matrix ( p. 10 ) with three angles and three phases 



that enters the interactions of neutrinos with the SM particles. In the MSSM, however, 
the phases ipA can be eliminated (by appropriate rotations of the superpartner fields) 
only if the slepton mass matrices and charged lepton mass matrices are simultaneously 
diagonal. In the general case, the phases ipA appear in the neutrino-sneutrino-neutralino 
and neutrino-chargino-charged slepton vertices. 



^In a general electroweak basis in which the leptonic Yukawa coupling is not necessarily diagonal the 
matrix U is given by the product U = e\v l where the rotations e E^e and // V lv diagonalize 
the mass matrices of charged leptons and neutrinos, respectively. 
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It is sometimes convenient to work with complex neutrino masses m,y^ in eq. ( p. 131 ). 
The phases ai^2 from eq. ( |2.10| ) are then absorbed into the masses: m^^ ^ = \m^^ ^ \ exp(— 2iai ^ 
In particular, when CP is conserved the phases ai^2 = 0, ±7r/2 reflecting different CP par- 
ities of different neutrinos, can be absorbed into the mass eigenvalues and make them real 
positive or negative. The MNS matrix is then real. Since only relative signs of the masses 
play a role, we will for definiteness fix our convention so that m^^g is always positive. Such 
a parametrization simplifies greatly the qualitative analysis of the renormalization group 
equations discussed in the next section. 

The value of the mass M is not known. One can expect that it is determined by a 
beyond the SM theory that provides the physical cut-off for the SM. The idea of Grand 
Unification and of the big desert would suggest very high value of M. However there are 
also other models (large extra dimensions) in which M is 0{1 TeV). In this review we 
assume that whatever the theory of neutrino masses is, the scale M is high enough to 
justify the description of the neutrino mass effects by the effective theory (SM or MSSM) 
with the single operator (|3.11|) added. 

Finally, one should also mention that in the MSSM neutrino mass can originate from 



R-parity violating interactions Such models are usually based on low energy mech- 
anisms for neutrino mass generation and will not be discussed in this review. 



4 Quantum corrections from the renormalization group 
evolution 

The neutrino mass parameters given at some scale M in the effective theory (SM or 
MSSM) supplemented by the operator (|3.11|) , determine the measurable quantities after 



inclusion of quantum corrections. The first category of corrections are those described by 
the evolution of the effective theory parameters from the scale M down to some low energy 
scale close to the electroweak scale. They depend on arbitrary powers of the large loga- 
rithms of the ratio M/Mz and are resummed to all orders of the perturbation expansion 
by means of the renormalization group equations (RGE). The second class of corrections 
are the so-called low energy threshold corrections. In the SM they connect the mass pa- 
rameters of the SM Lagrangian renormalized at the scale Mz to the fermion masses in the 
the effective theory obtained after the electroweak symmetry breaking (whose renormal- 
izable part is QCD+QED) renormalized at the same scale Mz- In the MSSM, threshold 
corrections include also the superpartner contributions due to their mass splittings and/or 
flavour violation in slepton mass matrices. The relative magnitude of the two types of cor- 
rections will be discussed later on. In this section we assume that the RG corrections are 
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the dominant ones and, after deriving the relevant RG equations, we discuss the potential 
role played by these corrections for the neutrino masses and mixing. 



4.1 RG equations for the CKM matrix 

Before we derive the RGE for the neutrino masses and mixing we discuss the technically 
easier case of the evolution of the Yukawa matrices and the CKM matrix. 

The renormalization group equations for the Yukawa couplings of the SM and the 
MSSM at one and two loop p3| , are well known. For convenience we reproduce 



the one loop RGE in Appendix B. They have a matrix structure which encodes the 
evolution of their eigenvalues and of all entries of the CKM mixing matrix. From the full 
RGE in the matrix form it is possible p5| , p^ ] to derive the RGE for the eigenvalues of 



the Yukawa matrices and for the CKM matrix. It is often very instructive to discuss the 
evolution of those Lagrangian parameters that can be directly determined from the data. 
For instance, the experimentally known information can be unambiguously extrapolated 
to any high scale. Moreover, certain qualitative features of the Yukawa coupling pattern 
necessary to reproduce the experimental data can be easier to understand. The derivation 
goes as follows. One writes down auxiliary RGEs for the matrices C/l,_R; D^^Ft defined in 
eq. (|3.4|) in the form: 

jD, = D,.i, ^n„-D,.l (4.1) 

where the matrices 5^'^ are antihermitean, in order to preserve unitarity of the matrices 
Dl,r, and Ul,r, and t = (l/lGvr^) ln((5/M^). One then looks for the evolution of the 
Ul^r and Dl,r matrices such that Yukawa matrices remain diagonal during the evolution. 
Differentiating the four relations (obtained from eqs. (|3.5|) ) 

diag(|/^, yl, y^) = U{YiY JJ ^ = U^^Y^YlUn 

diag(|/2, y'^, yl) = D{y\Y,Dl = D^Y^Y^Dn (4.2) 

and requiring the derivatives of the diagonal matrices on the rhs to be also diagonal 
matrices one gets the matrices s^^: 

2 2 

{^l)JI - o _,,2 1^^ ^ ydK^ \^l)JJ - U 

iJuj Uui K 

{e1)KL = -rf„444 E ylV'''*V'^, ie?)KK = (4.3) 

ydK ydL J 
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and similarly the matrices e^^\ 

tluj Uui K 

{sDkl = E yiy*V'', {e^nUK = (4.4) 

where Ud = du = —3/2 for the SM and Ud = du = 1 for the MSSM. The evolution of the 
CKM mixing matrix V is then given by the simple matrix equation 

^V = -elV + Vel (4.5) 

Because of the hierarchical pattern of quark masses the expressions ( [4.3|) for ej^^ can be 
simplified to 

{sDkl = -rf^y^V^^V^^ for K>L and (ef )xl = -{^^lk for K<L 
{e1)ji = -UdV'^V'^*yl for J > / and {e1)ji = -{el)u for J < J. (4.6) 

Taking next into account the hierarchy of the entries of the CKM matrix it is easy to find 
that, to a good approximation, the 2x2 submatrix describing the mixing of the first two 
generations as well as the element V^^ do not evolve, while the evolution of the remaining 
entries is universal and is given by 



V{t) = \/(0) exp |- [udvlit') + d^y^.it')] j (4.7) 

where V = V^^, V^^, V^^ or V^^. It follows that to one loop the Jarlskog invariant 
J = iYa(y'^'^V'^^V^^*V'^'^*) does not change. By using e.g. the standard parametrization 
of the CKM matrix it is also easy to derive the RGE for the three mixing angles and the 
CP violating phase.0 For the SM and the MSSM this has been done up to two-loops in 
ref. m\. 



4.2 RGEs for neutrino masses and mixing 

The Majorana mass term for neutrinos arises from the dimension 5 operator ( p.ll| ). The 
RGE for C^^ in the SM has been correctly computed only recently in ref. |2^ (previous 

^In the one loop approximation the running of the CKM matrix elements looks particularly simple 
in the Wolfenstein parametrization |28| (see e.g. (l6|) in which only the parameter A changes with the 
scale. 
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calculations 0, ^ have errors) and in the MSSM in refs. 0, 



dt 



-KC 



K 



{YlY:f C + C[y\Y, 



(4i 



where in the SM k = -3/2 and K = -Sg^ + 2Tr {syIY^ + ^Y^Y^ + YIY^) + 2A; the 



normalization of A is fixed by the Higgs self interaction: C, 



self = -^{H^Hf. In the MSSM, 
K = +1 and if = -Qgl - 2g'^ + 2Tr [^Y^Y^ . Since our initial conditions will be always 
given at the scale M, we have written eq. ( |4.8| ) top-down, i.e. with IGvr^t = ln(M/Q). 
The equation (^4.8|) is valid in any fiavour basis, in particular in the basis in which the 
Yukawa matrix Y^ is diagonal.0 In that basis, eq. ( [4.8|) simplifies to 



d 



dt 



C 



AB 



-KC 



AB 



K 



,,2 r^AB 



+ c 



AB„.2 



where y^^ are the eigenvalues of the hermitean matrix Y\Y ^ 
elegantly solved 



C{t) = IkJC{Q)J 



(4.9) 

In this form it can be 
(4.10) 



where J = diag(J, 



/,,, Ir) and 



Ik = exp ( - K{t')dt' 

ft 



(4.11) 



Note that any zero in the initial matrix C^^ is preserved by the RG evolution and that 
the phases of the C^^ entries do not evolve |3^. We also note that, since IrIc ~ 



IS 



the experimental constraint ( p.l4| ), with ml^ = (f ^/4M)C^"^(0) and = v'^{Mz) 
renormalized negligibly. 

Although the solution (^4.10|) to the RG equation for C^^ is simple, qualitative features 
of the running of the neutrino mass eigenvalues m,^^ and of the MNS matrix are often 
masked by the diagonalization procedure. In that approach it has to be performed after 

^°The MSSM RGEs given in ref. |Q allow to treat also the case in which squarks and/or gluino are 
much heavier than sleptons, charginos and neutralinos so that the decoupling procedure can be 
employed; there are then four different operators which mix with each other below the squark/gluino 
threshold. Above it one has (in the notation of ref. |3^) c'f' 
equations of ref. ||3l[| merge into the one quoted here. 



n ab 
^'-12 



n ab 
ZC21 



.ab 



C"" and the four 



llr 



This is because the matrix El defined in eq. (3.4) does not evolve; the matrices ef analogous to 



the ones defined in eq. (4.1) are identically zero because below the scale M, at which the operator (3.11) 



is generated, the RGE for the leptonic Yukawa coupling ( |B.iyB.5| ) do not depend on the couplings Y^, 
the only one that could change their matrix structure during the evolution. 
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the evolution. It is therefore useful to derive the RGE directly for the mass eigenvalues 
and mixing angles. This is done by using a trick similar as for the CKM matrix. One 
defines the auxiliary antihermitean matrix e'^ by the equation 



dt 



(4.12) 



where U satisfies eq. ( p.l3| ). Differentiating the equality C°'6°'^ = {U'^CU)"'^ and requiring 
that its rhs remains a diagonal matrix gives us the matrix e'^ . The only difference is that 
the resulting equation depends on {e^Y instead of {e^Y . Hence, one obtains two separate 
equations for Re(e'^) and Im(e'^): 



\ A 



(4.13) 



where 



- m„ - m 



(4.14) 



and of course Re(5'')'^'^ = Im(£^)"" = 0. For conserved CP {U real), those equations have 
been first derived in ref . . For a general complex matrix U they have been given in . 
The formulae for dU /dt written there is valid in a general electroweak basis; it reduces 
to ([4.13|) after setting P = Y\Y e = -P^ and passing to the basis in which P is diagonal. 
Thus, 



d 
dt 

d_ 

'dt 



U 



Aa 



Ah 



A,aReij2u'''*ylU 



Ba 



\ A 



\u 



B 



Aa\2 



+ Im 

Aba 



Ba 



B 



(4.15) 
(4.16) 



Eq. ( [4.15|) gives directly the running of the angles 6*12, 023, 6 



an 



if A parameterizing the matrix U 



Aa 



piifATjAa 
^ '^MNS 



13 and the phases 5, 
37| . In general, even if the phases 



are zero at some scale they will be generated during the evolution. However, it is easy 
to see that the differential equations for the parameters of the MNS matrix ( [2.1(J| ) do not 
depend on the ipA- Indeed, from eq. ( |4.13| ) it follows that does not depend on ipA] 



furthermore, since U 



Aa 



Aa 
MNS 



+ i(pAUyms\^ all factors e 



cancel out in eq. ( [4.12|) 



If two of the three eigenvalues, say m,j^ and niy^, are equal at some scale t there is the 
freedom in choosing the matrix f/(t), corresponding to the redefinition U{t) — > U{t) = 
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U{t)R where R is a rotation in the ab plane. For the evolution, R has to be fixed by the 
condition 

MEu''''*ylu^']it) = o (4.17) 

A 

SO that e'^ in nonsingular. This is particularly important for considering mass patterns 
with exact degeneracy at the scale M. The mixing matrix U is then ambiguous at that 
scale in the tree level approximation and becomes determined by quantum corrections, 
no matter how small. We shall illustrate this point in the next section. Note also that 
for conserved CP and neutrinos z/^ and z/f, of opposite CP parities such an ambiguity of 
rotations R does not exist. In such a case it is convenient to work with a real matrix U^"" 
and with the neutrino masses of opposite signs (see the comments at the end of Section 
3). The RG equation ( [4.15| ) is then nonsingular because Aat, instead of being divergent, 
vanishes. 

Eqs. ( [4.15| ) and (|4.16| ) are very convenient for a qualitative discussion of the impact 
of the RG evolution directly on the neutrino masses and mixing angles. As we shall see, 
several physical effects are in this approach more transparent than in the approach based 
on evolving the matrix C^^. 

4.3 Evolution of the neutrino masses 

The effects of the RG quantum corrections on the neutrino mass eigenvalues are simple 
and, except for a few special cases, not very interesting. The solution to eq. ( [4.16| ) reads 

m.„(t) = J^exp j\l{t')\U^^{t')\'dt'^ m,„(0). (4.18) 

Neglecting the small i/e and Yukawa couplings^ 

<it) = mliO)llexp (-4«: j\lit')\Usa\'it')dt'') (4.19) 

Factors Ik and k are different for the SM and the MSSM {k = —3/2 and +1, respectively). 
We see that the possibility of some change in the mass pattern caused by the evolution 
resides solely in the differences in the mixing matrix elements U^a and their RG running. 
In the MSSM where ^ (tan/^/lOO)^, tz ^ 0.12 for M = 10^° GeV and f/|„ typically 
varies between and 1/4 (except for U^^), the exponent is at most of order of e = 
y^log{M/Mz)/167i^ ^ tan2/3 x 10'^ < 2.5 x 10"^ for tan/3 < 50. In the SM, ~ iq-^ 

^^In the SM the Yukawa coupUngs are unambiguously determined at the electroweak scale by the 
charged lepton masses; in the MSSM (ye^)MSSM = (ye^)sM/ cos/3 where tan/3 = Vu/vd can vary from a 
few up to ~ 50. 
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and, consequently, e ~ 10 ^. We can then estimate the changes in the mass squared 
differences: 

Aml,{t) ^ mist) - ml^t) = Am2,(0) - (r^.m^JO) - Vbmlmt (4.20) 

where we have neglected Ik which is always close to 1, and the factors ria[b) > are 
typically in the range — 2, depending on the values of U^a factors and their evolution. 
Taking (for definiteness) Am^^(O) = 0, we see that the evolution of Am^^(t) is limited by 
m^^(0)e or m^^(0)e, i.e. by the value of the larger mass. 

We see that the RG evolution cannot change the pattern of masses. It may still be of 
importance for precision tests of various models, particularly those with partly degenerate 
or degenerate patterns. We recall here that such patterns have (at least approximately) 
degenerate two and three neutrinos, respectively, at the scale M. Matching the exper- 
imental Am^ and AM^ needs then some fine-tuning between the initial values and the 
RG quantum corrections. That point can, however, be meaningfully discussed only for 
each concrete model. 

The special cases are those with equal masses, m^^ = m^^ 7^ m^^ or m^^ = ml^ = ml^ 
at the scale M. They are called in this review two-fold and three-fold degeneracies, 
respectively, to be distinguishable from the more general partly degenerate or degenerate 
patterns. The interesting question we shall discuss in Subsection 4.5 and Section 5 is 
whether only quantum corrections can then explain the observed mass squared differences. 

In the next sections we shall discuss the RG quantum corrections to the neutrino 
mixing. We shall focus on interesting qualitative effects that depend only on the broad 
classification of the neutrino mass pattern, and do not depend on such details as whether 
the masses are equal or only approximately equal. 



4.4 Mixing of two neutrinos 

We first examine the mixing of two neutrinos, which was investigated in many papers 



3^ , p8t |39| , [40| , [4l|] . Strictly speaking, it could be physically relevant only if the atmospheric 
neutrino anomaly was due to the Un-Ur oscillations whereas the solar neutrino deficit 
resulted from the t'g-'^steriie oscillations (at present strongly disfavoured by the SNO and 
Superkamiokande data). Moreover the two 2x2 neutrino systems would have to be 
completely independent due to some particular texture of the U matrix. Nevertheless 
we will see, that the evolution of the 623 angle in the 3x3 scenario is, in some cases 
very similar to the evolution of the mixing angle of two neutrinos only. It is therefore 
instructive to discuss the 2x2 evolution and to compare it later with the more realistic 
3x3 mixing. 
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We begin with a real matrix C^^ and hence, a real U^'^: 

U^'^=h (4.21) 

where = sini), = cost?. We shall consider the same or different CP parities of the 
two neutrinos m^-^m^^ > or m^-^m^^ < 0, respectively. From eq. ( [4.15|) one finds the 
following equations for the mixing angle 

— = KA2i{y^^ - yeJsi)Ci) — = -KA2i[y^^ - Ve^s^^c^, (4.22) 



where n = —3/2 and +1 in the SM and the MSSM, respectively. Eqs. ( ^.22] ) give 

4 sin^ 2d = 2kA2i (yl - yt \ sin^ 2^ cos 2d (4.23) 
at V 2 1/ 

known from the literature^ ^ ^ Eq. (|]^) has a trivial fixed point (FP) 

at sin^ 2'(9 = but, contrary to the statements made in some papers, the maximal mixing 
€\r?2'd = 1 is not its FP: Although, naively, sin^ 2'i9(t) = 1 solves eq. ( [4.23| ), it is easy to 
see that s^{t) = ±l/\/2 does not solve eqs. ( |4.22|) .n 

It is clear from eq. ( |4.22| ) that, for fixed values of the Yukawa couplings y^^ and for 
some fixed evolution "time" tz = (l/lGvr^) ln(M/M2), the evolution of the mixing angle 
depends on the factor A21 = {my^+my^Y / Am\^. It is always small when the two neutrinos 
have opposite CP-parities i.e. when my^rriy^ < (as follows from eq. ( |4.18| ), m^^ cannot 
change sign during the evolution) or if the neutrino mass spectrum is hierarchical (m^^ <^ 



rriy or ml -C ml ). The evolution of the mixing angle is then negligible [B^, R^, H^, 



For m^^m^^ > the evolution of d can be significant, particularly if the neutrino masses 
are nearly degenerate so that \A2ie\ ^ 1 where e = (t/g^ ~ VeJ'^z- It can be checked that 
for |A2i|e > 3 the FP at sin^ 2'(? = is reached at the electroweak scale. 

In view of the maximal z/^-z/^ mixing needed to explain the atmospheric neutrino data, 
some attention was paid to the possibility of increasing d by the RG corrections from a 
small value at the M scale to a (nearly) maximal (l-i^l ~ 7r/4) at the electroweak scale 
3^ , |38| , |39| , As we said earlier, the maximal mixing is not a FP of eqs. ( [1.22| ), so 



such an "explanation" of large mixing can merely be due to a coincidence of the running 
"time" and the initial values of the angles and the masses. The value of |A2i|e must be 



^■^Its rhs is more frequently written in the equivalent form 2/^ ^^^ j^pii {Uez ~ Vei) ^i^^ 2i9cos^ 2-&. 
^■^For example, for constant 77 = nA^iiyl^ — Uei): eq- ( 4.22 ) is solved by s'^{t) = Sq/(sq + Cq^) where 



^ = exp(— 2ryi) and sq is the initial value of s^. It is then straightforward to check that sin^ 2?9(t) = 
sin^ 2??oC/(so + CqC)^ solves eq. (4.23) for any initial value sin^ 2i?o. Thus, for sin^ 2i?o = 1, eq. ( |42^ ) has 



two solutions: sin^ 2d{t) = 1 and the one given here but only the latter satisfies the underlying eq. ( 4.22|) 
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in the range such that the evolution is non-neghgible, but not strong enough to reach the 
FP. This can be most easily seen if we return to the solution ( |4.1CI| ) written in the form 



where hiit) = -^62 /-^ei = exp (^—k JoiVe^ — y1^){t')dt'^- One then has 

sin^ 2'dit) = (4.25) 

It is obvious that sm^ 2'dit) ^ 1 is obtained whenever C{t) ^ C^^(t). Since C'^'^it) 
evolves differently from C{t) (as ^ y1^) it is relatively easy to device the situation 
in which, at the initial scale M, sin^ 2'i9o is small (this requires \Cl' — Cg^| ^ 2|Co^|) 
and the evolution is such that at some lower scale C{t) = C'^'^{t) holds (this obviously 
requires Cq'^Cq' > 0). With the judicious choice of the C'^^(O) matrix elements and 



e = (yl^ — ye^)tz it is possible to obtain sin 2'(9 ^ 1 at the electroweak scale g^, |41 
Expressing Cg^ and Cq^ in terms of the neutrino masses and mixing angle at the initial 



scale M, the relevant condition Cq = ll^(tz) reads 

("^S^4o + <^4o) = + m(°)4j/|,(t^) (4.26) 

where s^^ = sin7?(0) and m^^J oc C"(0). It is clear that for ^ I'nT'^^l and s^g ~ 

(or ^ 0) satisfying the condition (f4.26|) requires rriy^m^^ > 0. On the other hand, 
for m^^rriy^ < the product Cq^Cq^ can be positive only if m^^ ^ "^^i "^^2 "^^i; 
which leads to a hierarchy Cq' > Cq"^ or Cq' < (7^2. Getting C2^(t) = C^^(t) requires 
then I21 ^ 1 or ^ 1. Thus, for m,j^my,^ < the evolution parameter e must be large, 
too large to be accommodated in realistic theories. This is in agreement with our earlier 
observation that for opposite CP parities of the two neutrinos and/or their hierarchical 
masses the RG evolution is very weak. 

Of course, since sin^ 219 = 1 is not the FP of the RGE, sin^ 2'd = 1 can hold only at 
one particular scale. 

Similar strategy can be also applied to analyze the general complex C^^ . Once the 
parameters of the U"^^ matrix 



''^ e-J(.: c jio e-«) (*-^^> 

are expressed explicitly in terms of the C^^ entries, the solution ( [4.1CI| ) allows to obtain 
analytic formula for sin^ 29{t) and to study its behaviour. The relevant formulae have been 
given in ref. (see also [0]). As could be expected, for fixed initial values of sin^ 26 and 
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the neutrino masses Im^,^ | and Imj^j 1; sin^ 29 obtained by the RG evolution interpolates 
smoothly between its value obtained for 02 = (i.e. m^-^m^^ > 0) and 02 = ±7i"/2 (i.e. 
with m^^mu2 < 0) 



4.5 Mixing of three neutrinos and fixed points 

For conserved CP, it is straightforward to derive from ([4.15|) the equations for the three 
independent mixing parameters S12, S23 and S13 0. Neglecting i/e and Yukawa 
couplings we get, both in the SM and in the MSSM: 

-512 = ~Ci2(Si2S23 — Ci2C23Si3)(— C12S23 — Si2C23Si3)/tA2i?/^ 

- Sl2Cl2C23Sl3(Sl2S23 - C12C23S13) nAsiy^ (4.28) 
+ Ci2C23Sl3(-Ci2S23 - Si2C23Sl3)KA32yr 

■S23 = ■Si2C23(Si2S23 — C12C23S13) K.A^iy'^ 

- Ci2C23(-Ci2S23 - Si2C23Si3)/t^32l/r (4-29) 

Sl3 = -Cl2C23Ci3(si2S23 - Ci2C23Si3)kA3i?/^ 

- Si2C23Ci3(-Ci2S23 - Sl2C23Sl3)ft;^32yr- (4.30) 

The evolution of the mixing angles can be classified into several universal types of be- 
haviour, depending on the magnitude of the factors Aab in eqs. (|4.28|44.30D . We note that. 



neglecting the small effects of mass evolution, all possible mass configurations in patterns 
listed at the end of Section 2 (except for the case of neutrinos (approximately) 
degenerate in mass and all having the same CP parity - to be discussed later) give one of 
the following four structures: 

a) A31 ^ A32 and IA31I ^ \A2i\ ^ 1 

b) A31 ^ A32 and \A2i\ > IA31I, \A2i\ > 1 
C) A32 ^ A21 ^ 0, 1^31 1 > 1 

d) A31 ^ A21 ^ 0, IA32I > 1 

For hierarchical masses and for partly degenerate structure with opposite CP parities 
of the (almost) degenerate neutrinos all Aab are 0{1). For partly degenerate pattern 
with same CP parities or for degenerate pattern, at most one of them is large. In the 
first case, it follows from eqs. ( [4.28|) -( p:.30|) and the value of e ranging from 10~^ in the 
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SM to 2.5 X 10 ^ in the MSSM with tan/? ^ 50 that the evolution of the angles is very 
weak |3^, 0. On the other hand, for one of the Aab factors sufficiently big so 



that \Aab€\ ^ 1, the angles evolve to an infrared quasi- fixed point (FP) [0, It is 
clear from eqs. (|4.28|) - ([4.30|) and from the parametrization ( 2.10 ) of the MNS matrix U, 



that depending on which Aab is large, the fixed points are either at f/31 = or U32 = 0. 
Before discussing the approach to those fixed points in more detail, we can already now 
summarize several qualitative conclusions. 

It is interesting to notice that in both fixed points we get the same relation between 
the mixing angles 

2^12 = .2% , J% with .^3 = - 1 ± Vl - sin^ 2^23 . (4.31) 

Thus, contrary to the 2x2 case, quantum corrections can give now an interesting, non- 
trivial fixed point relation. It is particularly interesting in the context of the present 
experimental indications for small 613 angle from CHOOZ and maximal atmospheric and 
solar mixing (sin^ 2^23 ~ 1 and sin^ 2^12 ~ 1). The relation ( [4.31| ) is inconsistent with such 



a pattern of mixing. We stress that quantum corrections summarized in RGBs eqs. (|4.28 ) 



( [4.3U| ), if large, always give ( 4.31| ). Thus, if the presently most likely pattern of mixing is 



confirmed experimentally, all mass patterns generating large quantum corrections through 
RGE are ruled out! 

Special cases easy to consider (still before a detailed study of the approach to the fixed 
points) are exact degeneracies at the scale M: ml = ml 7^ ml or ml = ml = ml . As 



follows from the discussion surrounding eq. ([4.17|) , for the same CP parities, m^-^ = m^^ 



my.^ the angles must satisfy the FP relation ( [4.31| ) already at 
the scale M. If this is ruled out by experiment, then the two-fold degeneracy needs 
my^ = — m,^2 and the three- fold degeneracy is ruled out. This conclusion holds (both in 
the SM and MSSM) under the assumption that the RGE corrections are the dominant 
ones (see next section for other possibilities). For m^^ = —m^^ and m^^ ^ ^^35 with 
\ml^ — m^U ~ AM^, we can ask if quantum corrections can explain Am^ = \ml^ — m^J. 
The answer to this question is positive in the MSSM (in the SM see [^) and will 



be discussed in more detail in Section 5. Another point worth a discussion is what if 
experiment will eventually be consistent with the FP relation (^4.31|) . Two-fold degeneracy 
is easy: we need m^^ = m^^, \ml^ —^l^l ~ AM^ and (as discussed in Section 5) quantum 
corrections can explain Am^ . With three-fold degeneracy, we necessarily have at least one 
pair of neutrinos with the same CP parities but the question is whether we can explain 
both Am^ and AM^ by the discussed here class of quantum corrections. The answer is 
negative |]51| , ^ This can be explained as follows: from eq. (|4.19|) we would need 
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|f/3i| ~ |f^32| (to keep ml_^ ^ '^^j) during the entire evolution and, at the same time, 
l^^ssl ~ \U32\ ~ C^(l) to generate sufficiently large AM^. With 613 ~ the ffist condition 
implies sin^ 26'i2 ~ 1 and the second sin^ 26'23 ~ 1 (If^ssl ~ |f^32| can be at most ~ 1/2). 
This is however incompatible with the relation ( 4.31| ) which should be satisfied! So, even 
if experiment was consistent with ( [4.31 ), the three-fold degenerate mass spectrum would 
still be unacceptable since AM^ and Am^ could not be explained by quantum corrections 
(always under the assumption about the dominance of the RG corrections). 

Let us now discuss the approach to the fixed point and concentrate first on the struc- 
ture b) which can be realized for partly degenerate or degenerate patterns, with same CP 
parities of z/i and 1^2. Since A31 ^ A32 and \A2i \ ^ l^sil, equations (|4.28|-|4.30|) reduce to: 



S12 
S23 

Sl3 



-Ci2(Si2S23 - Ci2C23Sl3)(-Ci2S23 - S12C23S13) ^21?/^ - C12S23C23S13 A32l/r ) 
S23C23^32l/^, 



(4.32) 



023513013^32?/^ 



(we take MSSM, with k = 1). The equation for S23 is the same as in the 2x2 scenario 
discussed in the previous subsection. The evolution of ^23 in the two cases is formally 
not identical as the evolution of the mass factor ^32 depends now also on the remaining 
mixing angles 612 and 613. Nevertheless, the qualitative behaviour of S23 is similar because 
in most cases the scale dependence of ^32 can be neglected. Denoting 



= exp 1^ 2A^2{t')yl{t')dt'^ ^ exp (-2A32(0)e) 
the solution for 5^3 (t) reads 



'23 



'23'' 



(0)/ 4(0) + 4(0)^. 



and yields 



sin^ 2^23(t) = ir sin^ 2^23(0)/ 5^3(0) + 4(0)^. 



The solution for 5^3 can also be given in a closed form: 



'13 



W = 4(0)/ {4(0) + 4(0) [4(0) + 4(o)e.] } 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



Thus, since |y432| ~ 1, the evolution of both, S23 and S13 is very weak |3^, 0. For 
example, the effect of the running for sin^ 2^atm is a 2.5% change for extreme value of 
tan 13 ^ 50. 

For |v42ie| ^ 1 (recall we consider ^^^2(1) ^ |Am^| and my^rriy^ > 0), the evolution is 
towards one of the two approximate fixed points of the RG equation for su. One can 
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easily check, for instance, by considering the equation for d tan 9 12 / dt, that for A21 > 
(i.e. for Am^ > 0) the point f/31 = is the UV fixed point and U32 = is the IR fixed 
point. For A21 < (i.e. for Am^ < 0) the situation is reversed. It is also interesting 
to notice that in the limit S13 = we can follow analytically the approach to the fixed 
points. In this approximation 

S12 = si2c\2slj,A2iyl (4.37) 
and the solution for is of the form (|4.34|) , with ^23(023) — > 512(012) and[^ 



- e = exp (- 2sl,{t')A2,{t')yl{t')dt'^ . (4.38) 

For y42i > 0, in the top-down running, the factor ^' ^ exponentially with decreasing 
the scale Q Mz i.e. for growing t oc A2iy'^\og[M / Q) and, consequently, we obtain 
Si2{t) = ±1 (depending on its initial sign) and approach IR fixed point at Us2 = 0. 
Changing the signs of the right hand sides of eqs. (|4.32| ) one can see that in the bottom- 
up evolution we approach Si2(t) ~ exponentially, i.e. the UV fixed point at f/31 = 0. 
For A21 < we get the reversed situation, in accord with our general expectations. 

It is interesting to estimate the values of A21 and tan/3, for which the approach to the 
IR fixed points is seen. In ref. we estimated that for approaching the fixed point during 
the evolution in the range (M, Mz) with M ^ 10^° GeV one needs A2ie{M) > 3, i.e. for 
tan/3 = 20 one needs rriu-^^ ~ mi,^ > 10~"^ eV for Am^ ~ 10~^° eV^ and m^^^ ^ m^^ > 0.01 
eV for Am^ > 10~® eV^. The qualitative change from a negligible evolution to the FP 
behaviour at the electroweak scale is abrupt and occurs in the small range 0.5 < \Aabi\ ~ 3. 

Finally we note that, from the point of view of the initial conditions at the scale 
M, the UV fixed point looks not realistic as the neglected muon Yukawa coupling 
quickly destabilizes it during the evolution. We conclude that for hierarchical and partly 
degenerate mass patterns the evolution of the mixing angles is either very mild or shows 
(for \A21\e ^ 3) a fixed point behaviour. 

The evolution of mixing angles in the degenerate case, m^^ ^ ml^ ~ m^^ ~ (9(AM^) 
or larger, partly falls into the same classes of behaviour. Indeed, as long as ^31 ^ ^32 
with 1^31 1 < 0(1) and IA21I ^ 1, the angles evolve according to the same equations 
One can easily identify the mass patterns of the degenerate case that fall into 



^'"'Eqs. (4.37) and ( [4.38[ ), after solving for S23(t) and in the approximation of constant A21 and A 
give: 



23, 



'12 



it) - 42(0)/ {sUo) + clM [clM + slM^T^''^^'^}- 
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this category: the necessary condition is that m^^ and m,^^ are of the same sign. For the 
evolutions of S12 we then closely follow the two possibilities, depending on the sign of A21, 
discussed for the partial degeneracy with 171^^^171^^ > 0. We simply note that larger values 
of \A2i\ are generic for the present case and the approach to the fixed points is faster. 
However, the evolution of S23 and S13 are guaranteed to be mild only if m,,-^ and Tfi are 
negative. For positive and my^i i-e. when all masses have the same CP parity, we can 
have |v432e| > 3 for tan/5 > 40 (due to the bound (|2.14|) \A^2A ~ 3 cannot be realized for 
tan/3 < 40). The angle ^23 behaves then as the angle d of the 2x2 mixing discussed in 
Section 4.3 and, according to the solution (|4.34|) , S23 is exponentially focused to the stable 
FP S23(t) = or S23(t) = ±1, depending on the sign of ^32, and on the direction of the 
evolution. The angle S13 behaves in a similar way except that, as follows from eq. ( [4.36| ), 
it does not reach the value 5^3 = 1 when S23 1 (due to the presence of the factor C23 in 
its RGE, the evolution of S13 is then "frozen"). We conclude that in the regime in which 
the approach to the fixed points is relevant, the pattern with approximately degenerate 
neutrino masses and the same all three CP parities is not acceptable. 

As we said earlier, the approach to the FP behaviour is abrupt as a function of 1^316 1 
1^3261 . However, in the small transition region of the values of ^31 and ^32, in agreement 
with our discussion in the previous subsection, it is possible to chose the initial condition 
for S23 so to get sin^ 2^23 = 1, ■S13 ~ at the electroweak scale. This was exploited in 
ref. |^5[ as a possible mean to obtain maximal atmospheric neutrino mixing from the 
initially small sin^26'23. Since |v42ie| is always much larger than |A32e|, the evolution of 
S12 is then such that the FP at f/32 = or f/31 = is quickly reached. Thus, a realistic 
solution, with maximal 6'23 and S13 ~ 0, has sin^ 2^12 ~ and, moreover, the scheme 



cannot work unless there is an extreme fine tuning of the initial parameters |37 



The remaining degenerate mass patterns can be classified according to the relations 
A32 pa A21 ^ or A31 ^ A21 ~ 0. Consider first A21 ~ A32 0, i.e. rriu^ ^ —rriy^ ~ rriy^. 
The equations for the evolution of the mixing angles can be approximated as 

512 = -Sl2Cl2C23Sl3(Sl2S23 - Cl2C23Sl3)A3lZ/r> 

S23 = Si2C2(Si2S23 - C12 023513)^31?/^, (4.39) 

513 = -Cl2C23C3(si2S23 - Ci2C23Si3)v43i?/^. 

These equations exhibit IR quasi-fixed point behaviour for v43ie ^ — 1, corresponding to 
t/31 = 0. As before, at the fixed point the angles satisfy the relation S13 = tan^i2tan6'23. 
Since S12 is proportional to S12 and suppressed by S13, the running of S12 is weak. The 
IR fixed point is reached due to strong running of S23 and S13. For y43ie ^ 1, S23(M^) is 
strongly focused at ±1. Thus, the mass and tan/3 configurations leading to A2,ie ^ 1 are 
unacceptable. For ^31 ^ A21 ~ and ^32 < we get IR fixed point in U32 = 0. 
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The RG running of the mixing of three neutrinos for a complex matrix U, i.e. non- 
zero phase 6 and ai, 02 7^ modulo 7r/2 with all masses m^^ positive by definition, has 



been also investigated in the literature For initially degenerate Ua and Ub the mixing 
pattern is determined by the condition Re (u^"-U^^*^ = Re (f^MNS^MNs) = playing in 
the general complex case the role of the FP condition. Qualitatively, our conclusions 
remain unchanged also in this case. For example, starting with m^^ = wi^^, sin^ 2^23 = 
S13 ~ and arbitrary phases ai, 0^2 at the scale M, the RG evolution leads to the relation 

m 

sin^ 2eu{Mz) = sin^ 2^i2(M) sin^ai - 02) + 0{sl^) (4.40) 

We see therefore that maximal solar mixing can again be obtained only if the change of 
612 during the running is negiligible, i.e. for ai — 02 ~ 7r/2 at the scale M. 

The running of the masses is always given by eq. ( |4.18|) and is similar as for a real 



matrix U. The conclusion that the triple degeneracy of neutrinos at the scale M, |mj,J = 
Irriiy^l = cannot lead to an acceptable pattern of masses and mixing also remains 

valid. 

In summary, with all \Aab€\ ^ 0.5 the evolution of the mixing is negligible. For the 
mass configurations such that at least one \Aab\ ^ 1 and \Aab€\ ^ 3 the infrared fixed 
points are reached during the evolution, independently of further details of the mass 
matrices. However, only for \A2ie\ > 3 and IA31I, IA32I ~ 1, or for ^31(^32)6 ^ —3 and 
^31(^32) ~ A21 ~ the evolution is consistent with a large atmospheric mixing angle at 
low energy. The mass configurations m^-^ ^ m^,^ ^ rrii,^ and ±mj,^ ^nii,^ niy.^ with 
^31(^32) > also lead to the infrared fixed points but at the same time the atmospheric 
mixing angle converges to zero. We also note that for |v42ie| > 3 only S12 runs to assure 
the fixed point relation, so the initial values for S23 and S13 have to be close to their 
experimental values already at the scale M. For ^316 < —3 or ^326 < —3, S23 and S13 
evolve strongly and the evolution of S12 is weak. The IR fixed point relation ( |4.31| ) is 



always one equation for three angles. Insisting on a large atmospheric mixing angle, it 
correlates small (as follows from CHOOZ) ^13 angle with a small solar mixing angle. The 
fixed point solution makes the low energy angles dependent on only two, instead of in 
general three, initial conditions at the scale M. We conclude that quantum corrections 
encoded in the RG running of mixing angles may have dramatic impact on their physical 
values if the mass pattern is partial degeneracy or degeneracy. If the bimaximal mixing 
solution was confirmed, all the mass patterns leading to the FP would be ruled out, unless 
the low energy threshold corrections change the results. 
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5 Low energy threshold corrections 



5.1 Threshold corrections in the SM 



The Wilson coefficient C^^ of the dimension 5 AL = 2 operator ( |3.11| ) is a renormahzed 
parameter of the effective theory Lagrangian. Integrating its RGEs from the high scale 
M down to some scale Q Mz allows to resum potentially large corrections involving 
\n{M/Q) to all orders of the perturbation expansion. However, since the low energy scale 
Q is not a priori determined by any physical requirement (apart from the condition Q ^ 
Mz), the neutrino masses and mixing angles computed in the tree-level approximation 
from the Wilson coefficient C^^{Q) do depend (albeit weakly) on the actual choice of Q. 
This dependence can be removed by computing masses and mixing angles in the one-loop 
approximation in the MS scheme and with Q as the renormalization scale. 

In general, to compute the physical neutrino masses and mixing one has ffist to perform 
the RG evolution of the entire matrix U and neutrino masses from the scale M down to the 
scale Mz, include threshold corrections and subsequently rediagonalize the resulting mass 
matrix ml^^^°°^ by an additional unitary matrix U'. The physical matrix ?7mns (whose 



rAa 



elements are probed in neutrino experiments) is then given as e^'^'^U^^g 
where U is the matrix obtained from the RG evolution. 

In the basis in which neutrino masses are diagonal, we have 



{U ■ U') 



i\Aa 



m 



1— loop 
ab 



(5.1) 



where 



A,B 



u 



Bb 



(5.2) 



In the SM one finds 156 



7-th 
^AB 



167r2 2 M^ 



+ ... 



(5.3) 



/M^. Note that since g'^ml /2M\ 



ul^, the coefficient of \n.{l/Q) 



where 

agrees with the coefficient k of the y^^C^^ + C^'^yl^ term in the RGE 



for the 



SM. This confirms the correctness of the recent re-derivation |^ of the SM RGE. The 
remaining corrections indicated by dots in ( p.3| ) affect only the overall scale of the neutrino 
masses and therefore are not interesting in view of the unspecified magnitude of the mass 
M in eq. (^1^ . 

As could be expected, in the SM the nontrivial part (|5.3| ) of the low energy threshold 
corrections that changes the structure of the mass matrix, and hence of the matrix U"^^ 



IS 
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small and always negligible compared to the quantum effects described by the RG running, 
which are enhanced by large logarithm of the ratio M/Mz- The correction ( p.3| ) can be 
most easily taken into account by stopping the RG evolution of the Wilson coefficient 
C^^ at the scale Q = Mwe'^^'^'^. 



5.2 Threshold corrections in the MSSM 

In the supersymmetric model, apart from stabilizing the results obtained from the RG 
analysis, the low energy threshold corrections can be as or more important than the RG 
evolution and may have very important physical consequences for the neutrino masses and 



mixing angles Before we discuss some physical examples, it is worthwhile to 

adapt our calculational procedure to the new situation, so that the RG evolution effects 
and the threshold corrections can be treated on equal footing. Although the general 
procedure outlined earlier is in principle correct irrespectively of the relative magnitude 
of the RG and threshold effects, in practice it can mask simple qualitative features if the 
rediagonalization due to threshold corrections is not a small perturbation. 

We observe that, to a very good accuracy, in the solution (|4.1CI| ) to the RGE the factors 



Je^ ( [4.11| ) can be approximated by 



h.^l-nf" yljt' ^ 1 - fl (5.4) 





With this approximation (which can fail only for unrealistically large values of tan (3) all 
quantum effects of the physics below the scale M, the RG running as well as the low 
energy threshold corrections, can be described by a single formula (|5.1|) . The and 
the matrix U are now the neutrino mass eigenvalues and the neutrino mixing matrix, 
respectively, at the scale M, and the factors 1°'^ are given by eq. ( |5.2| )with 

1% - 1% - II ^ Iab- (5.5) 



Again, the case of a real matrix U is particularly easy because then the formula (|5.1|) 
can be written as 

^a^'°°" = m.J'^' + (m,„ + m.J/'^^ = m^J'^'' + Am,5 (5.6) 

with the right hand side symmetric and real i.e. Hermitean. One can then use the formal 
perturbation calculus of quantum mechanics (see e.g. |^) to find corrections to neutrino 
masses and mixing angles. Of course, if the threshold corrections are absent]^ one has 
Iab ~ 5^^Ia with |J^| > |/^| > 

^^If the threshold corrections are universal, i.e. if 1^^^ = I^^d"^^ ^ they can be absorbed into the overall 
scale of neutrino masses and do not influence neither the mixing angles nor the ratio AM^/Ato^. 
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Several physically interesting situations can be discussed. In the following we shall 
mainly focus on the scenarios with equal three or at least two masses at the scale M. 
Although the discussion of the mixing angles does not depend on whether the masses 
are equal or only approximately equal, the former case is much more constrained for the 
masses themselves and therefore more interesting. 

5.2.1 Three-fold degeneracy and flavour diagonal corrections 

In the MSSM the complete formulae for 1^^^ are lengthy due to additional contributions of 
if^ and chargino/charged slepton and neutralino/sneutrino loops but we do not need 
them for our discussion. First of all, it is interesting to investigate whether it is possible 
to generate the right mass squared differences starting from degenerate neutrinos at the 
scale M: |mj,J = \my^ \ = Irriu^l (we again allow for different CP parities of the neutrinos). 
From the arguments presented in Section 4 it follows that this is impossible if the factors 
Iab are dominated by the RG effects, i.e. when Iab ~ S'^^Ia with ^ |/^|, |/e|. The 
degeneracy Ansatz is however potentially interesting as it might help to understand the 
bimaximal mixing. Also, the neutrino masses of order of few eV are necessary if the 
neutrinos are to play a role of the hot component of the Dark Matter.[] It is therefore 
interesting to see if some other corrections can split the mass squares. 

We start with the flavour diagonal threshold corrections: Iab = ^^^^a and consider 
first the mass pattern m^^ = —m^j^ = = rriu- The formula (|5.6|) now takes the form 



/1 + 2UIJA 2UaiUa3Ia\ 



1— loop 



(5.7) 



-I-2UX2IA 

\2UaiUa3Ia 1 + 2UI^IaJ 

The corrections to neutrino masses can be calculated perturbatively. Because rriui = "it^us, 
one apphes here the perturbation calculus to the case with degeneracy of the unperturbed 
"energy levels" [0. Instead of solving the "secular" equation for the corrected eigenval- 
ues, it is however better to exploit the freedom of an arbitrary rotation in the plane (13): 
U —>■ UR13. This freedom can be used to diagonalize the perturbation by requiring 

T.UaiUa3Ia = (5.8) 
a 

which fixes the matrix U after taking into account the perturbation. The correction Amab 

to the zeroth order neutrino mass matrix becomes then diagonal and the neutrino masses 

^^The hot component of the Dark Matter was previously needed to account for formation of largest scale 
structures in the Universe. At present, in view of the observational evidence for a significant contributions 
of the cosmological constant (or another form of dark energy as e.g. quintessence) to the energy density 
of the Universe, the hot component of the Dark Matter seems no longer necessary (but is not excluded). 



29 



are given by 

\m,J = m,{l + 2UljA). (5.9) 

Let us now assume [^] that Jg ^ 0, = 1^ = 0. The condition (|5.8| ) reduces then to 
UuUis = 0. The experimentally viable solution is Uis = 0, i.e. S13 = 0. Note that this 



relation replaces the FP relation ( 4.31 ) of Section 4. For the neutrino masses one finds 
Am2i = -4m^ cos^ 26*12/6, Amga = -4m^Si2/e- (5.10) 

It is therefore possible to generate AM^ ^ Am^ provided the solar mixing is nearly 
maximal. The same result is obtained starting from — mj,^ = m,^^ = rriy^ = my. Thus, if 
the underlying theory gives one of these two degeneracies together with the bimaximal 
mixing pattern, the low energy threshold corrections can be responsible for assuring S13 = 
and generating the correct mass squared differences. 

Of course it is unrealistic to expect only one correction I a to be nonzero. Note however, 
that the results for the mixing angles and the ratio I\M'^ j lS.vp? are not altered by shifting 
all corrections I a by an overall additive constant: Ia = ~ I ■ Indeed, using 



unitarity of the matrix ?7, the off-diagonal entries in the rhs of eq. (5^) can be always 
written as X^a UAaUAt^A (o- 7^ b) with arbitrary I. For the diagonal entries instead one 
has 1 + 2 Ea UUa = 1 + 2/ + 2 Ea UIJa ~ (1 + 2/)(l + 2 Ea UIJa)- Hence, up to 
higher order terms, eq. ( p.6| ) can be rewritten with = J2a UAaUAh^A and rriy replaced 
by fhy = (1 + 2I)my. Thus, = Jr 7^ is equivalent to = Jr = 0. 

It follows that |/e| ^ |/^|, equivalent to |/e| ^ |/t| 7^ 0, = 0, should give only 
small correction to the result ( |5.10| ) and S13 = 0. For rriy^ = —rriy^ = m^^, solving eq. (|5.8| ) 
(with /e = le, 1^1 = and 1^ = Ir) one finds |Q that 

C12 

Aml^ ^ -Amlhsl^ (5.11) 



Am2i ^ -Amlle 



cos 2^12(1 - sl^r) + (1 + 2c?2)si3 



where r = Ir/Ie- Thus, obtaining small but non-zero angle ^13 is also possible. To obtain 

the experimentally favoured value of AM^ one must have \m1Ie\ ~ 1.6 x 10"'^ eV^ i.e. 

Ve ~ ^ 10^^ and ~ Ir for rri;^ in the eV range. In the MSSM the required hierarchy 

|/e| ^ |/t-|. It ~ Ifi can arise from the low-energy threshold corrections only as a result 

of non-universality of the left-handed charged slepton masses (and sneutrino masses) .p^ 

-"^^In the loop approximation to I*^^ it was estimated in ref. that in order to get /g ~ lO^'^ and 
positive one needs A/g^ « 1.7Af^^ (and -M^^^ « Mf^). However, from the fuU expression for /^'^ it follows 
p6[ that /e ~ 10~^ and positive can be achieved only for relatively light first chargino (toci ~ 300 GeV) 
and Mji^ w (1.2 — 1.6)Mg^. For Mg^ « (1.2 — 1.6)M^^, which can be realized in the inverted hierarchy 
models |5p, one gets le ~ —10^"^ i.e. positive Am|]^ for (almost) maximal solar mixing. 
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Moreover, for exactly maximal solar mixing, obtaining Am? appropriate for the LAMSW 
solution requires s^g ~ 10"^^"^-*. This is consistent with the CHOOZ data and is realized 
for r < 0.1, i.e. |Jt- — /^| ^ 10~^. For tan/? < 3 such a small difference between the 
corrections I-r and can be due to the RG effects in /^^ and does not require any mass 
splitting between /i/, and tl. For larger values of tan/3, for which I^^ — IJf > 10~^, some 
conspiracy between — due to M^^ < M^^ and /^^ — IJf is required. The amount 
of the necessary fine tuning grows, of course, with tan/5. For the VO solution to the 
solar neutrino problem, with sin^ 26*12 = 1, a more severe fine tuning is always necessary: 
Am^ ~ O{10~^^ eV^) requires | J^— /^| < 10^'' which is always smaller than I^^—P^ > 10~^. 
Thus obtaining the correct Aim? for this solution requires a cancellation at least to one 
part per hundred between the RG and threshold effects. For An? ~ (9(10~^° eV^) one 
gets ^ 10~^. Similar results are obtained also for — m,^^ = rrii,^ = m^,^. 



It can be checked [60| that, for the initially degenerate neutrino masses, the patterns 
m^^ = —niiy^ = m^.^ or —m^-^ = m^^ = m^^ and |Je| ^ \Ir\, It ~ are the only ones 
that can produce the required mass squared splittings with flavour diagonal threshold 
corrections. The necessary condition is that the underlying theory valid above the scale 



M gives the bimaximal mixing in the basis in which eq. ( |5.8| ) is satisfied [^H 

5.2.2 Three-fold degeneracy and flavour non-diagonal corrections 

Qualitatively new possibilities for mixing angles and for splitting initially degenerate 
neutrinos originate from nonzero off-diagonal elements of the corrections in eq. (|5.2|) 
[61 , 3C]. This is possible in the MSSM if the slepton mass matrices are not diagonal in 



the flavour space in the basis in which the leptonic Yukawa couplings are diagonal. The 
amount of flavour mixing is then best quantified in terms of the so-called mass insertions 
||6^ defined as the ratio of the off-diagonal (in flavour space) elements of the charged 
sleptons mass squared matrices to some average charged slepton mass squared.|3 Current 
limits on such mass insertions following from the non-observation of the decays /i ^ 67 
etc. are not very stringent: Only the insertion mixing left(right)-handed fl with 
right (left)-handed e has to be smaller than O{10~^). The insertions 6}^j^ and 6]^^ causing 
similar, chirality changing, e ^ f and fl ^ f transitions are bounded by ~ 0.5 and ^ 0.1, 
respectively, for slepton masses ~ 500 GeV. Bounds on the chirality preserving insertions 
"^LL ' ^RR exist only for the e ^ ft transition and are ~ 0.2 for ~ 500 GeV. The other 
chirality preserving mass insertions are practically unrestricted. 

^^Due to the underlying S'[/l(2) symmetry, mass insertions inducing transitions between left-handed 
sleptons of different generations are always accompanied by the insertions inducing similar transitions 
between sneutrinos. 
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If the mass insertions are non-vanishing, there is a flavour non-diagonal contribution 
to I^^. In general it takes the form 



^AB^TT—, E S^^hxY{M~,^^,mc^,mr,^) (5.12) 

-■-"^ X,Y=L,R 



where Hll, Hrr and hiR are some functions of the chargino and neutralino masses ^«Cj , 
rrijq. and of some average mass M^^^ of charged sleptons. The largest values of /^*^ can 
be obtained for relatively light charginos, heavy sleptons (~ 1 TeV) and M2/ ^ ^ — 1 and 
reach /i^^l/IStt^ < few x lO^"' (^r_r and h^R are always smaller) For comparable 

chargino and slepton masses one has /ill/IGtt^ ^ 2 x 10~^. In principle the mass insertion 
approximation should fail for |5xyl ~ 0.1. In practice it works as an order of magnitude 
estimate even for \5xy\ ~ 1 (^^e error is then of order 25%). 

In the presence of non-zero mass insertions and for niy^ = m^^ (= —m^J the condition 
for vanishing of the appropriate off-diagonal entry of the correction to the zeroth order 
mass matrix reads 

J2UAaUBblAB = 0. (5.13) 

A,B 

Consider first the situation in which the single correction I^^, J*^ or dominates over 



all other corrections. The condition (^.13|) gives then relations between the mixing angles 



that are different from (|4.31|) and are listed in Table |l| Only three of the nine 

possibilities are compatible with the bimaximal mixing: dominant for m^-^ 



or = ''^!^3 and dominant for m^-^ = m^^. For initially degenerate neutrinos, the 
latter combination gives wrong relation Am^ ^ 2AM^. (Other six combinations leading 
through ( ^.13| ) to sin^ 26i2 ~ sin^ 2^13 also give bad relations Am^ ^ AM^ or AM^ = 
0.) The former two are however interesting giving Am^ ^ — 4m^ cos 26'i2 sin 2^23-^^t ^^^1 
AM^ ^ 4ml{l + CI2) sin 26*23/*^. Obtaining AM^ ^ 3.2 x 10"^ eV^ is therefore possible 
but only for > 1 eV and > 0.5 i.e. for rather large flavour mixing in the slepton 
mass matrices. 

Similarly as for a non-zero Jg correction, the solar mass squared difference Am^ <^ 
AM^ can be generated either by an appropriately tuned departure of the angle |6'i2| from 
maximal value tc/A or by another, hierarchically smaller, correction / 7^ 0. It has been 
demonstrated , that including on the top of the dominant correction a hierarchically 



smaller perturbation in the form of either J^, P {I^ does not work) or I^^, I^^ allows to 
split mj^j = ~m^2 even for exactly bimaximal mixing 



20 



We note |60|, that with a large le diagonal perturbation discussed in the Subsection 5.2.1, Am <^ 

-th 



AM^ (and 7^ s^g <C 1) can be also induced by or I^^ instead of / = (small 7*'^ 7^ would 
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Table 1: Relations of the FP type between the mixing angles for dominant correction 



m 



m 



f2 



m 



m 



f3 



— 72 — 

1^ 



1^ 



7-th 



Sl3 = cot 26^12 cot 6[ 



23 



cot 6^12 tan 6*23 



tan 6*12 tan 6 



23 



rth 



Si3 = - cot 2^12 tan 6123 



l-2s 



cot 612 cot ^23 



l-2s 



tan 612 cot ^23 



rth 



2si3 

l + s2 



tan 26*12 tan 2^23 



S13 = — tan 6*12 cot 26; 
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S13 = cot 6*12 cot 29\ 
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In all those cases of two hierarchically different corrections the important difference 
between the two alternatives: perturbation by a diagonal correction I a or perturbation 
by an off-diagonal correction J^*^ is that, for exactly bimaximal mixing, in the former 
case Am? oc r^, whereas in the latter Aim? oc r only (si3 ~ r in all cases), where r <^ 1 
is the ratio of the smaller to larger correction as in ( |5.11 ) |30|. Therefore, similarly as 
for |/e| ^ obtaining Am? appropriate for the VO solution with dominance and 
hierarchically smaller ler or Jg^ would require some tuning of slepton masses to cancel too 
large a contribution of /^^ to /,-. 

Finally, we remind the reader that the relations listed in Table 1 remain approximately 
valid when the equalities of the masses are relaxed and replaced by the corresponding 
approximate degeneracies. The listed relations play then the role of the FP relations 
discussed in Section 4. They are satisfied at the electroweak scale irrespectively of the 
initial values of the angles. The observed Am^ and AM^ can be obtained by adjusting the 
initial values of only approximately equal masses. The role of the threshold corrections 
Pab is then the same as the role of the RG corrections. If large enough, they give one of 
the "fixed points" relations of Table 1. 



5.2.3 Two-fold degeneracy and threshold corrections 



We can also discuss the effect of threshold corrections in the case of the two-fold degeneracy 

the Ami 9 off- 



2 2 
m,, = m 



ml > ml^ or m^^ 



For m 



-m 



V2 



diagonal entry in eq. (|5.6| ) automatically vanishes but the correction matrix Amah as 
the whole needs not be diagonal. In the first order of the perturbation calculus p7| , 
the neutrino masses are then given by eq. ( ^.9|) and receive also further corrections of 
order O ((Amafe)^/max(|mi,3|, |m,y|)). The mixing angles also receive corrections of order 



-rs23(c23) {r = Ilt.jie) and Am^ 



not split TOjyj^ and m^^ for sin 29i2 — 1). This leads to S13 
4m^ sin20i2C23(s23)^e)i(r)- Obtaining Ato^ ~ 0(10"'* eV^) requires, for m,, 
^1^(t) ^ 10~(''~^^ which is possible in the MSSM and gives r ~ 10~^^^^^ i.e. acceptable S13. 



1 eV the correction 
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O (Ama6/max(|mj,3 1, |m,^|)), i.e. small if the hierarchy of neutrino masses is large. If only 
flavour diagonal threshold corrections are present the solar mass squared splitting is given 
by0 



Am2i ~ 



4m. 



(5.14) 



The interesting aspect of this situation is that the (generically) dominant first order 
contribution to Am^ proportional to Ir = I^^ + Z**^ vanishes for f/|. 



f/|i, i.e. for si3 



±(523/023) (ci2 ± Si2)/(ci2 T ■512), a special case of which is the bimaximal mixing solution 
with Si3 = (also the second term in eq. ( |5.14| ) then vanishes). This has been discussed 
in refs. ^ p6l in connection with the possibility of realizing the VO solution to the 
solar neutrino problem, within the inversely hierarchical pattern m^^ <^ ml ~ 3 x 10~^ 
eV^. Note that for generic mixing angles the RG corrections would always give too large 
a Am^. 

Since the second order correction (oc /^) to Am^ is proportional to m^^g, the right Am^ 
for the VO solution can be obtained by appropriately tuning m^^ > and/or due to the 
threshold effects. 



As discussed in [061, with m 



-m,. 



3 X 10 ^ eV^, the correct Am^ 



for the LAMSW (or SAMSW) solutions can be obtained from the RG running. From 
eq. ( |5.14j) we get 



Am ^ 4m, 



S23 cos 2^12 + S13 sin 2^23 sin 29 



12 



(5.15) 



and for the LAMSW solution, with |si3| ~ 0.1, Am^ ~ (9(10 ^) eV^ can be obtained for 
tan/3 > 30 (tan/? > 10 for the SAMSW solution). For smaller values of tan/3 one has to 



investigate potentially larger threshold corrections, which give |60| 



Am^ ^ Ami 



Am^ 



4m 



sin2^12(c23/e;-S23/, 



th^ 



For Am^ ~ (9(10 ^) eV^ one needs therefore corrections / 

rth 



th 



LAMSW 
SAMSW. 

10-3 - 
rth 



(5.16) 



definitely too big a value for If^ in the MSSM and also impossible for I^^ and J*^ (at least 



10 ^ which is 

th 
e(/^)r| 



within the validity of the mass insertion approximation, it would require | > 1, 

Our discussion here does not cover fully the potential role of quantum corrections for 
the two- fold degeneracy pattern. Several other possibilities do exist, depending on the 
chosen solar solution and we refer the reader to the literature |^1], |5^ for further 
details. 



^We choose to work with = — 
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6 Conclusions 



For every theory of neutrino masses, for a meaningful comparison with experimental 
information, it is necessary to discuss quantum corrections. This is relatively easy if the 
effective low energy theory is the SM or its supersymmetric extension and the neutrino 
masses enter via the effective operator ( |3.11 ). In Sections 4 and 5 we reviewed the 
formalism for including quantum correction in those cases. Next we applied that formalism 
to several different classes of the neutrino mass sector at the scale M, hypothetically 
given by the theory of neutrino masses. They split into two broad groups. The first one is 
characterized by small quantum corrections that may eventually be important for precision 
tests of the future theory. However, they are irrelevant at the level of present, qualitative 
considerations. Here belongs the hierarchical pattern and the partly degenerate pattern 
with opposite CP parities of the (almost) degenerate neutrinos. For those mass matrices, 
quantum corrections cannot substantially alter the structures present at the scale M, so 
the agreement with experimental data has to be assured by the boundary values at M. 
Equal masses = —mu2, at the scale M are possible. Quantum corrections can explain 
the observed Am^. 

The other group consists of partly degenerate pattern with the same CP parities 
of the (almost) degenerate neutrinos and of the (approximately) degenerate structures. 
Large quantum corrections can originate either from the RG evolution or from low energy 
threshold corrections (in the MSSM). They never change qualitatively the mass eigenvalue 
pattern, although they may explain their observed splitting. However, large quantum 
corrections always lead to a "fixed point" relation for the mixing angles. The sufficient 
condition is that for at least one pair of neutrinos {m^^ + m^J/{m^^ — m^^) ^1. It is 
interesting that the transition from small (qualitatively irrelevant) quantum corrections 
to the "fixed point" behaviour is very abrupt. So, to a good approximation there are 
those two and only two physical situations. 

A fixed point relation is always one equation for three angles (if CP is conserved) and 
makes their low energy values dependent on only two boundary conditions for the angles 
at the scale M. If the dominant quantum corrections come from the RG evolution (i.e. 
originate from the large r Yukawa coupling) or also from several configurations of the low 
energy threshold corrections, the "fixed point" relation links small ^13 angle (constrained 
by the CHOOZ experiment) to a small angle responsible for the mixing of solar neutrinos. 
If presently favoured bimaximal mixing was confirmed by future experimental data, all 
mass patterns leading to large RG corrections and the regions of the MSSM parameter 
space leading to the same "fixed point" relation would be ruled out. 

There are, however, other "fixed point" relations, generated by some special sfermion 
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mass configurations and/or by flavour non- conserving effects in the slepton sector, that are 
consistent with bimaximal mixing and small ^13 angle. Such solutions are phenomenologi- 
cally interesting as, at the same time, they explain the observed mass squared differences, 
as the effect of quantum corrections, with degenerate spectrum at the scale M.0 However, 
the simple Ansatz at the scale M needs a deeper theoretical justification. 

Quantum corrections do not explain the origin of the neutrino masses and do not 
replace its theory. Nevertheless, they are important piece of the overall picture. They 
will constrain strongly the acceptable mass structures once the experimental ambiguities 
are resolved. 
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Appendix A 

In Section 3, instead of the familiar Dirac bispinors we used the two-component Weyl 
spinor notation. It is particularily convenient for dealing with Majorana particles. Here 
we explain this notation shortly. More details can be found in modern textbooks on QFT 
(see e.g. |l^, |64 



In four dimensions, the Lorentz group or more precisely its covering group SL{2,C) 
has two non-equivalent complex two-dimensional representations denoted as (1/2,0) and 
(0,1/2). The Grassmann fields (or fermionic field operators) transforming according to 
these representations are conventionally written as and and called^ left- and right- 
handed spinors, respectively. Since the complex conjugation of a left-handed spinor Aq, 

^^It has also been discussed in the hterature ^ that quantum corrections generated in the fuU 
theory above the scale M, can break the degeneracy aheady at the scale M and lead to correct mass 
squared differences and mixing angles at the electroweak scale. For an interesting link between neutrino 



masses and flavour changing processes see ref. 63 . 

^■^Somewhat incorrectly; properly they should be called left- and right-chiral. 



36 



transforms as a right-handed one {{Xa)* ~ A^), the fermion content of any Lagrangian 
can be specified by hsting only the left-handed spinors used for its construction. 

If the two left-handed fields Aq, and X/3 transform as representations R and R*, respec- 
tively under the final unbroken symmetry group (global or local) of the theory, they can 
be combined to form a Dirac bispinor: 

^(A) = Q;), v^(A) = (x",A^), (A.l) 

transforming as R and R*, respectively. The raising and lowering of Weyl spinor indices 
is done with the help of the antisymmetric tensors e"^, e^/?, e^^ and e"'^: 

A- = e°%, A„ = AV> Xa = e^pXp, r = Xpe^''- (A.2) 
The two kinetic terms for A and x can be then rewritten in the familiar form 

^kin = iXaf'd^X + ix^^d^x = iXa^d^X + ixcr'^d^x + (total der) = iip^x)l^d^ip(x) (A.3) 
where the Dirac matrices in the Weyl representation are constructed as 

with cr^ = {I,(t), d"^' = {I,—cr) (cr's are the Pauli matrices). For such a pair of Weyl 
fields also a Dirac mass term can be constructed 

Cmass = -m (A"Xa + "^AaX") = -'^^{A)^(A) (A. 5) 

If the field A (x) has no left-handed partner transforming in the complex conjugate rep- 
resentation R* (R), it is convenient to introduce chiral Dirac bispinors 

For chiral Dirac bispinors e.g. V'(a)l = Pl4'{x)l etc., where Pl = {I — 7^)/2. Note that 
(see eq.O)> 

V^=(0 A^), ^=(x" 0). (A.7) 

The typical Yukawa coupling of a scalar field (p in the representation R^ and two left- 
handed Weyl spinors A and x transforming as representations R\ and R^, respectively 
(such that 1 G R(i,x RxX R^) can be written as (we omit the Clebsch-Gordan coefficients) 

-CYuk = ~Y<j)Xx - Y*<j)^Xx 

= -Y(t)i){^^)Ri){^x)L - Y*(f)'^ij^x)L^ix)R (A.8) 
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The Yukawa part of the SM Lagrangian ( |3.1| ) and the SM mass terms ( ^.31 ) are the example 
of ( [A .81) and ( |A.5| ), respectively, with fields u^, df^ and e'^ playing the role of x, and q and 
I (or u, d and e) playing the role of A. 

Finally, Weyl spinor fields which are singlets of all unbroken symmetries of the 
theory can form 4-component Majorana bispinors 



a 



^(A)Maj= Q^j. (A.9) 

Of course in this case ipix) = C'V'p^) = V'(^a) which means that the field is self-conjugate. 
One also has (up to raising or lowering indices) 



A = ^/'(A)L = ^(A)i? and A = ^/'(A)i? = V^(a)l- (A.IO) 
For such a field a Majorana mass term can be formed 

-^Maj = -^m(AA + AA) = -^m (^(A)i?^(A)L + ^{A)LV'(A)i?) 

= -^mV^(A)V^(A) = -^mijf^^Cijf^x) (A.ll) 

where C is the charge conjugation matrix. The Majorana mass term ( |3.8| ) is precisely of 
this form with u'^ playing the role of A and the Yukawa coupling ( |3.7| ) is usually written 

as 

A>CYuk = -eijHi^(^^c)RY^ipQ)L - eijH*ip^x)LYlip(^c)ji (A. 12) 



Appendix B 



In this Appendix we recall the well known RGEs for Yukawa coupling matrices defined in 
eq. (13) in the SM and in the MSSM. In the SM they read p5|, |24|: 



where 



dt 



d 

dt 
d_ 

dt 



Y^ 
Yd 
Y. 



- \9l ~f/y + T+\ [YiY^ - Y\Y, 
-89l - U -^g'y + T + l {YlY, - YIY^ 



4 

9 2 15 

92 



4 



' ^^91 + T+IyIY^ 



IGtt^ 



(B.l) 



(B.2) 



T = Tr {SYiYu + SY^Y^ + YIY, 



(B.3) 
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and the gauge couplings 5^3, (72 and gy evolve according to 

d 



dt 



9i = kQi 



3,2,y 



(B.4) 



with 63 = -7, 62 = -19/6 and by = 41/6. 
In the MSSM one finds instead: 



dt' 



d 

dt 
d_ 

dt' 



Yd 



Yd 
Y. 



—91 - H - ^^gl + Tr (3n>^«) + 3n^« + Y\Y, 
-y^73 - ^gl - \gl + Tr (3^^!-, + Y\Y + 'iY\Y , + YIY ^ 
-^gl - 3^7y + Tr [?>Y\Ya + YIY + ?,YIY 



(B.5) 



and the factors hi change to 63 = —3, 62 = +1 and hy = 11- 

For completeness we give here also the RGEs above the scale M i.e. for the theory 
whose set of fermion fields includes additional three SUl{2) x f/y(l) singlet neutrino fields 
u'^^ (y4 = 1, 2, 3) and whose Lagrangian is identical to the one for the effective theory valid 
below M except for the Yukawa interaction (|3.7|) and the Majorana mass terms ( p.8|) . 

If the theory above the M scale extends the SM, then T given in eq. ( B.3 ) has to be 
replaced by 

(B.6) 



T = Tr (3YiYu + 31^^^, + YIY, + YIY, 



the last equation in (|B.1|) should be replaced by 



dt' 



9 

'4 
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g2 



-gy 



+ t + -(yIy,-yIy, 



and the neutrino Yukawa matrix RGE reads: 



^Y 

dt •■ 



■lgl-lg'y + T-l{YlY,-YlY, 
In addition, the Majorana mass matrix also runs |B2|: 



(B.7) 



(B.8) 



dt "^^'^ 



Maj- 



(B.9) 



If the low energy theory is the MSSM, then ||39|, ^ in the first equation of ( |B.5|) 



Tr i'iYiYu) ^ Tr (?,YiY^ + yIy , 



the last equation of 
d 

^Ye = Ye 

dt 



rsl ) is replaced by 

-3gl - 3gl^ + Tr (3^^^^ + YIY,) + 3YIY, + Y\Y , 
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(B.IO) 



[B.ll) 



and the neutrino Yukawa matrix RGE reads: 



dt 



?>gl -gl + l^T (srir.. + YlY,) + ^yIy, + YIY,] . (B.12) 



Finally, the Majorana mass matrix running is dictated by |^ 



(B.13) 
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